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Abstract 



A test for the Gaussian and exponential Ansatze for the nonperturbative parts of the coefficient 
7-h ■ functions, j) non P eTt - and x>^ onpert \ which parametrize the gauge-invariant bilocal correlator of the 
field strength tensors in the stochastic vacuum model of QCD, is proposed. It is based on the 
! evaluation of the heavy-quark condensate within this model by making use of the world-line 
formalism and equating the obtained result to the one following directly from the QCD Lagrangian. 
This yields a certain relation between p nonpert - (0) and , D 1 onpert '(0), which is further compared with 
an analogous relation between these quantities known from the existing lattice data. Such a 



comparison leads to the conclusion that at the distances smaller than the correlation length of the 
vacuum, Gaussian Ansatz is more suitable than the exponential one. 
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Nowadays, stochastic vacuum model (SVM) ^ [| is recognized to be one of the most 
powerful tools for the investigation of both perturbative and nonperturbative phenomena in QCD. 
Although the predictive power of this model is very high, its field-theoretical status requires 
further justifications. In fact, an exact derivation of field correlators, which are up to now used 
in this model as a dynamical input, from the QCD Lagrangian is extremely desirable. Once 
being performed, such a calculation would shed a new light on an interpolation between the 
perturbative and nonperturbative effects in QCD. That was the main reason for various authors 
to address this problem analytically both in QCD M and other models possessing the confining 
phase || . However, despite these efforts, the exact form of the coordinate dependence of the bilocal 
gauge-invariant correlator of field strength tensors in QCD (which plays the major role in SVM) 
remains unknown. The most reliable result known from evaluations of the bilocal correlator on 
the lattice [|, [7|, || (see Refs. || [ToJ] for a review of the recent progress in the lattice calculations 
of field strength correlators) is that the nonperturbative parts of the two coefficient functions 
parametrizing this correlator, £> non P ert - an d 2)° onpert ', rapidly decrease at the distance, which is 
usually referred to as the correlation length of the vacuum, T g . The latter one is equal to 0.13 fm 
for the S77(2)-case and is about 0.22 fm for the SU (3)-case 0. 

The two standard nonperturbative Ansatze for 2) nonpert - and x>^ onpeTt \ used both in lattice 
investigations and in phenomenological applications of SVM to the evaluation of various QCD 
processes, are the Gaussian and the exponential ones. In the present Letter, we propose a non- 
trivial analytical test of these Ansatze. It is based on the evaluation of the heavy-quark condensate 
by making use of the world-line formalism [|II|] (see Ref. for a recent review) and further com- 



parison of the obtained result with the one, which follows directly from the QCD Lagrangian [13 
Such a comparison then leads to some relations between V nojipeTt -(0) and T>™ npert (0) for both 
Ansatze. Among those, as we shall eventually see, the one obtained for the Gaussian Ansatz satis- 
fies the existing lattice data j7|, || |9|, [Hj concerning such a relation much better than the other one. 
Then, since the typical sizes of the heavy-quark trajectories in the problem under study are of the 
order of T g , this leads to the conclusion that at the distances smaller than T g , Gaussian Ansatz is 
more suitable than the exponential one. This result is important for future lattice simulations as 
well as for the applications of SVM to the high-energy hadron scattering 
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Our method of derivation of the quark condensate is based on the well known formula 
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where V is the four- volume of observation and m is the quark mass. Next, the average (. . .) Aa 
on the R.H.S. of Eq. (|]) is implied w.r.t. the gluodynamics action in the Euclidean space-time, 
| / d 4 x fe) 2 , with a = 1, . . . , iV c 2 - 1 and F® u = - d v A* + gf^A^A^ being the Yang-Mills 

is the averaged one-loop quark effective action 



field strength tensor. Finally in Eq. ([I]), \ T 
(i.e. one-loop quark self-energy) defined as [[II], [L 
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Here, the subscripts P and A denote the periodicity properties of the respective path-integrals, 
ip p s are antiperiodic Grassmann functions, and = A p T a with T a 's standing for the generators 
of the SU(N c )-group in the fundamental representation. We have also adopted the standard 
normalization conditions (r[0]) Aa = and {W(0)) Aa = 1, where 

(W(C)) A% = ^trPexp (-ig J drA,x)j ^ (3) 

is just the Wilson loop, which is defined at the contour C parametrized by the vector x M (r). It 
turns out that the gauge-field dependence of Eq. (Q) can be reduced to that of the Wilson loop @ 
only. That is because, as it has been demonstrated in Ref. [15|| , the spin part of the world-line 
action can be rewritten by means of the operator of the area derivative of the Wilson loop and 
becomes separated from the average (. . .) Aa \ 
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Within the SVM, the Wilson loop @ can further we rewritten by virtue of the non-Abelian 
Stokes theorem and the cumulant expansion as follows ||1], [| 
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Here ((00')) Aa = (00') Aa — (0) Aa (0') Aa , and E[C] is a certain surface bounded by the contour C 

fi fj, fj, fj, 

and parametrized by the vector z^) with £ = £ 2 ) standing for the 2D-coordinate. This surface 
is usually chosen to be the one of the minimal area for a given contour C . We have also denoted for 
brevity z = z(l;), z' = z(!;') and introduced the notation Q(z, z') for the parallel transporter factor 

taken along the straight line joining the points z' and z: Q(z, z') = V exp i^—ig J A^{u)du^}j . 

Next, in a derivation of Eq. (j^), the so-called bilocal approximation has been employed, ac- 
cording to which the irreducible gauge-invariant bilocal field strength correlator (cumulant) 
((F^ u (z)^(z, z')F\ p (z')§(z' , z))) Aa dominates over all the cumulants of higher orders (see e.g. 

Ref. |16|] for details of a derivation of Eq. (^|) and discussion of related approximations). Finally, it 
is worth commenting that the factor 1/2! on the R.H.S. of Eq. (§) is simply due to the cumulant 
expansion, whereas the factor 1/4 is due to the (non-Abelian) Stokes theorem with the usual 
agreement on the summation over all the indices (not only over those, among which the first one 
is smaller than the second, used in Refs. [[j], ^ ^ |lf|). 

Let us further parametrize the bilocal cumulant according to the SVM by the two 
renormalization-group invariant coefficient functions T> and T>\ as follows 

((F, v (z)<f>(z, z')F Xp (z'Mz', z))) A = l J^u\{8 lxX 5 vp - 5 w 8 vX ) V ((z - z'f) + 



3 



1 

+ 2 



d d 
— ((z - z)\8 vp -(z- z') p 5 uX ) + gj- ((z - z') p 5 pX ~{z- z)\5 PiP ) 



VtUz-z' 



(6) 



Here, 1n c xN c is the unity matrix, and we have used the standard normalization condition for 

the functions V and T>\ || with Af = 2 4(p(o)+Pi(o/)' • ^ virtue of the (usual) Stokes theorem, 
the contribution to the double surface integral standing in the argument of the exponent on the 
R.H.S. of Eq. (^|), which emerges from the Pi-part of the cumulant @, can be rewritten as a 
boundary term. After that, we eventually arrive at the following expression for the Wilson loop 
within the SVM 
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Next, one can extract the volume factor V from the heavy-quark self-energy (0) by splitting 



the coordinate x p (t) into the center-of-mass and the relative coordinate fl^| , x p (t) = x p + y^r), 



where the center-of-mass is defined as follows x p = ^ J drXnir). The empty integration over x 

o 

then obviously yields the factor V . After that, owing to the fact that 5x p (t) = dyJr), we can 



also replace 



5o>i/(2i(t)) 
8 2 



in Eq. (m by 



<5o>„(j/(t)) 



, which is possible due to the formula |L 
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All the above mentioned considerations lead to the following intermediate expression for the 
heavy-quark self-energy (0) 
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Here, the contour T is parametrized by the vector y p {r), and the surface S[r], parametrized by 
the vector u> M (£), is bounded by this contour. 

Next, it can straightforwardly be shown that in the heavy-quark limit under study the typical 
quark trajectories are small. Indeed, let us consider the free part of the bosonic sector of the 
world-line action standing in the exponent on the R.H.S. of Eq. (P), 
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Here, in the last equality we have performed a rescaling r = |, T old = Among all the 

reparametrization transformations t — > er(i), % > 0, it is convenient to choose the proper-time 
parametrization t = — . Here, s is the proper length of the contour T so that T — — , where 
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L = L[r] = J ds = J daJy^ia) is just the length of Y. Within this parametrization, J dty 2 (t) = 

r ^in v o 

m J ds ( ) = since ( ) = 1 by the definition of the proper time. Therefore Sf ree = 

mL, which means that the typical quark trajectories are so that L < — , i.e. they are really small 
in the heavy-quark limit. 

However, according to the general concepts of the SVM, the distances smaller than the corre- 
lation length of the vacuum, T g , are forbidden for a test quark, i.e. L cannot be arbitrarily small. 
That is because SVM is an effective low-energy theory of QCD with T~ x playing the role of the 
UV momentum cutoff. Therefore, infinitely small contours C (whose contribution to the Wilson 
loop (0) would obviously dominate if they are allowed) are forbidden within the SVM. Moreover, 
C's should be so that not only their lengths obey the inequality L > T g , but even among such 
contours those which lie inside the sphere Sj^^) are forbidden. We conclude that within the 
SVM, typical heavy-quark trajectories are located from outside to the circle of the radius T g /2 and 
with the center placed at x, nearby to it (in the sense that their lengths L's obey the inequality 
L < —). Therefore, the correlation length of the vacuum is related to the area S inside such a 
trajectory according to the formula 

T 2 g * ~S (9) 

with a good accuracy. We also see that for any two points y^r) and y^r'), the following inequality 
holds 

\y(r)-y(r')\<T g . (10) 

Note in particular that since typical sizes of Wilson loops in the problem under study are ap- 
proximately equal to T g , our final conclusion on the advantage of the Gaussian Ansatz w.r.t. the 
exponential one will be valid at the distances \x\ <T g . 

Let us now evaluate heavy-quark condensate (p with the Gaussian and exponential Ansatze 
for the nonperturbative parts 2) non P ert - and x>" onper ' of the functions T> and T>\. The first of these 

2 2 
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Ansatze reads T> (x 2 ) = T>(0)e Tfl , T>\ (x 2 ) = T>i(0)e T f , where from now on we shall denote for 
brevity V nonpcrt - by V and £>° on P ert - by V x . By making use of the result of Ref. [1181 for the string 



tension (see also Ref. ]16[ for a detailed derivation), we get due to Eq. (§) the following leading 



P-function contribution 

2 J J da^{w')V ({w - w') 2 ) ~ AT 2 V{0) J d 2 te~ t2 ■ S ~ 16V(0)S 2 . 
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For further purposes let us rewrite this result as 8T'(0)S 2 tl/ , where = S M ^[r] = S^fC] = 
§ y^dy v is the tensor area associated with the contour Y (or C). Here, we have used the fact that 
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for contours under study, which slightly deviate from the circle, S 2 = \E 2 . By virtue of Eq. (|TUD, 
we also have 



i r (y- y >)2 r r 



dy, f dy'^ / dt V 1 (t) = V 1 (0)T 2 & dy, f dy'^e 
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This leads to the following expression for the Wilson loop standing on the R.H.S. of Eq. 
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with a s = |-. Note that the behaviour of small Wilson loops of the type e const - s was for the 



first time obtained in Ref. [@]. 

In the same way one can treat the exponential Ansatz, D (x 2 ) 



£>(0)e T « , £>i (x 2 
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Here in the evaluation of the V- and "Di-dependent parts, Eqs. ([]) and flTOP, respectively, have 
been applied. We see that for the exponential Ansatz, Eq. (|TT|) remains valid with the replacement 

"^Gauss * "^exp 2vAQ auss . 

In order to proceed with the evaluation of the path- integral (Bl), it is useful to linearize the 
quadratic E M „-dependence of Eq. (JTTJ) by introducing the integration over an auxiliary antisym- 
metric tensor field. Clearly, since E^ v depends on the contour T as a whole, this field will be 
space-time independent. Introducing for brevity the common notation A for both ^.Qauss and 
Ajxp, we have 
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Note that only in this equation we have emphasized explicitly the summation over indices in order 
to avoid possible misleadings. In what follows, in the expressions of the type O^O' we shall as 
usual assume the summation over all the values of indices. 

The one-loop heavy-quark self-energy (^) has now reduced to that in the constant 
field B^, which should eventually be averaged over. Indeed, since due to the Stokes 



theorem, B^H^ = J da illl (w)B lll/ (w) } the following equality holds (cf. Ref. |[L5| ) 
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Eq. (|TTD into Eq. (|8|) the following expression 
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Taking now into account the known one-loop expression for the Euler-Heisenberg Lagrangian in 
spinor QED, i.e. one-loop electron effective action in a constant background field (see e.g. [12"]), 
we have 



Vy^ / exp 



2 ^ iivV iiVv 



(4vrT) 



(4vrT) 



T 2 abcot(aT) coth(6T) - 1 



Here, the standard notations were adopted: a 2 
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with E = i {B 1A , S 24 , £34), H = (^23, Sis, B l2 ). 



Next, due to the factor e~ m T in Eq. (12), in the heavy-quark limit under study only small 
values of the proper time T are sufficient, and the R.H.S. of Eq. ( |TB"| ) can be expanded in powers of 
T. Then, the leading term of this expansion, which in the expansion of the fermionic determinant 
corresponds to the diagram with two external legs of the .B^-field, reads T 2 abcot(aT) coth(6T) — 
1 = (b 2 — a 2 ) + O (T 4 (E ■ H) 2 ). This leads to the following expression for the heavy-quark 
self-energy (|T^) 
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Next, the pole at D = 4 emerging in the integral over proper time, J dTT 
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mD j) 4 r(2- can be subtracted by making use of the MS'-prescription (Clearly, the addi- 
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tional constant, which appears if we use e.g. the MS'-prescription, does not contribute to the 
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quark condensate, since the derivative of this constant w.r.t. m in Eq. ([!]) vanishes.). This finally 
yields 
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with A standing for the UV momentum cutoff. Owing to Eq. ([[]), at N c — 3 we arrive at the 
following value of the heavy-quark condensate (M^ip) = which should be equal to the result 

following directly from the QCD Lagrangian (i.e. from the corresponding triangle diagram) Hl3| 
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For A = ^4causs, this yields the following relation between the 



A" 



nonperturbative parts of the functions V and T>i at the origin: Pi(0) = |£>(0). Similarly, for 
A = A cxp , we obtain Pi(0) = 5T>(0). These relations should be compared with the results of 
those lattice measurements 0, ||, [|, |10[| , whose x 2 is minimal. Such results can be summarized by 
the table, following from the data reviewed in Ref. [TO] , which has the form 



Theory 


Pi(0)/P(0) 


Quenched approximation at the distances 


0.22 ±0.05 


0.1 fm < \x\ < 1 fm 




Quenched approximation at the distances 


0.29 ±0.13 


0.4fm < \x \ < 1 fm 




The same as in the previous line, but with 


0.34 ±0.04 


the perturbative-like part of the fit fixed to 




zero 




Full QCD with the quark mass (in lattice 


0.13 ±0.08 


units) equal to 0.01 




The same as in the previous line, but with 


0.13 ±0.07 


the quark mass equal to 0.02 




Average over the above values 


0.20 ±0.10 



This table indicates that the result obtained within the Gaussian Ansatz is closer to the existing 
lattice data than the other one. 

In conclusion of the present Letter, we have proposed a nontrivial analytical test for the 
Gaussian and exponential Ansatze for the nonperturbative parts of two coefficient functions, which 
parametrize the gauge-invariant bilocal field strength correlator in QCD within the SVM. It was 
based on the evaluation of a heavy-quark condensate within this model by making use of the 
world-line formalism and further comparison of the obtained result with the exact one, following 
directly from the QCD Lagrangian. The outcome of this investigation suggests that the Gaussian 
Ansatz satisfies the existing lattice data better than the exponential one. This means that at the 
distances |x| < T g , the former Ansatz is more favourable than the latter one. Such a conclusion is 
important for further applications of the SVM to the evaluation of various QCD processes as well 
as for the lattice measurements of field correlators. 
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